function drawvec=drawntrunc(m,n,lims); 
%function drawvec=normtrrnd(m,n,lims);
% Lims  Limits 
% Results is an m x n Matrix of draws of 
% x ~ N (0,sigma)  |x| < lims
% % Modified 6/11/02
drawvec=zeros(m,n); 
num=1; 
while num <=m; 
   accept=0;
   while accept<1; 
      draw=randn(1,n);
      if abs(draw) < lims; 
         accept=1; 
      else;
         accept=0; 
      end;       
   end;
   drawvec(num,:)=draw;
   clear draw
   num=num+1; 
end;